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This work tests the statistical mechanical theory of hydrophobic interactions, isolates consequences of excluded 
volume interactions, and obtains B 2 for those purposes. Cavity methods that are particularly appropriate for 
study of hydrophobic interactions between atomic-size hard spheres in liquid water are developed and applied 
to test aspects of the Pratt- Chandler (PC) theory that have not been tested. Contact hydrophobic interactions 
between Ar-size hard-spheres in water are significantly more attractive than predicted by the PC theory. The 
corresponding results for the osmotic second virial coefficient are attractive (B 2 < 0), and more attractive 
with increasing temperature (AU2/AT < 0) in the temperature range 300 K< T <360 K. This information 
has not been available previously, but is essential for development of the molecular-scale statistical mechanical 
theory of hydrophobic interactions, particularly for better definition of the role of attractive intermolecular 
interactions associated with the solutes. 



I. INTRODUCTION 

Hydrophobic interactions are universally acknowl- 
edged as fundamental contributions to the stability of 
folded or aggregated biomolecular structures in water. 
But hydrophobic interactions are also expected to be- 
come more favorable with increasing temperature for 
physiological temperatures. On this basis, hydrophobic 
interactions are principally entropic interactions. 

Hydrophobic interactions can then be described as fa- 
vorable for aggregation and endothermic at moderate 
temperatures. The osmotic second virial coefficient 




= lim B 2 [R] (1) 

it— >oo 

is the solution thermodynamic metric for assessment of 
attractive and repulsive character of hydrophobic inter- 
actions between an AA solute pair. Here t/AA (f) is the 
usual radial distribution function of AA pairs at infinite 
dilution. Positive values of B 2 raise the osmotic pressure, 
and indicate preponderance of repulsive effects. Negative 
values of B 2 lower the osmotic pressure and characterize 
interactions that are attractive on balance. 

The entropic aspect of hydrophobic interactions re- 
quires statistical thermodynamics for explanation. Since 
the molecular theory of hydrophobic interactions thus re- 
quires specified intermolecular interactions and defensi- 
ble statistical mechanics, the theory of hydrophobic inter- 
actions has been intermittent and only partially success- 
ful. The approximate Pratt-Chandler (PC) theorjff' was 
the first prediction of molecular-scale <?aa (?") for inert 
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atom solutes in water and thus the first molecular-scale 
prediction of B 2 . Though the PC theory predicted at- 
tractive and repulsive features in unprecedented detail, 
it did not straightforwardly conform to the expectation 
that hydrophobic interactions as expressed by B 2 should 
be attractive and endothermic. The PC theory was im- 
mediately contro ver sial!^ Experiments for benzene, and 
perfluorobenzene^l disagreed with the PC theory for 
atomic-size hard-sphere solutes. Explanations for the 
discrepancy were suggested^] but the underlying con- 
troversy has persisted. 

One challenge for addressing this controversy is that 
the integrated quantity B 2 , and particularly the R — > 00 
limit, is a subtle target for molecular simulation calcu- 
lations. Another challenge is that this problem requires 
analysis of temperature (T) dependences in a limited T 
range. 

But further challenges remain, and a summary of 
the substantial efforts to resolve this issue was given 
recently.^ The broad conclusion is that simulation tests 
of the theory have differed enough in details that a con- 
clusive test of the PC theory for hydrophobic interac- 
tions has been elusive. For example, B 2 s for the specific 
cases of Ar and CH4 solutes have been estimated^ to be 
remarkably small, due evidently to substantial cancella- 
tion between repulsive and attractive force effects. In the 
same context,^ -B2 for Kr(aq) has been estimated to be 
repulsive (positive). Clearly, attractive and repulsive in- 
teractions can play conflicting roles with the consequence 
that merely realistic simulation of a case of interest, e.g., 
CH4, might not provide helpful physical conclusions. In- 
deed, the case of atomic-size hard-sphere solutes has not 
been treated specifically because hard-sphere models are 
inconvenient within common molecular dynamics tools. 

As another example, solute non-pairwise interac- 
tions associated with polarizabilities were of simula- 
tion interest for a timeP^ Those interests have waned 
inconclusivelyP^S] undoubtedly because of a lack of a 
molecular theory that could sharpen the questions and 
make the answers more permanent. Such differences are 
indeed broadly expected to be details, but in this context 
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FIG. 1. A configuration of 5x10 water molecules (red and 
silver) together with the inclusions (blue spheres) identified 
by nt = 2x10 trial placements of a hard sphere with dis- 
tance of closest approach to an oxygen atom of 0.31 nm. This 
size corresponds approximately to an Ar solute for which the 
van der Waals length parameter cta is about 0.34 nmP thus 
adopting 0.31 nm - 0.17 nm = 0.14 nm as a van der Waals 
contact radius of the water oxygen atom. Hard sphere so- 
lutes of this size are about maximal for water-oxygen contact 
densityP 



those details have been features of controversy. 

In parallel with simulation efforts, the foundation of 
the molecular theory of hydrophobic effects has under- 
gone a surprizing renovation! 10 * 13 ! The renovated the- 
ory directly exploits molecular simulation data and sub- 
stantially amends the PC theory, but gives unexpected 
support for some of the conceptual ingredients of that 
theory.^ For example, the concern that the PC theory 
neglected molecular orientational structure of liquid wa- 
ter is now generally recognized as not well founded. The 
same can be said about the scale d-par ticle approaches 
applicable to the present problem! 10 ! 15 ! The renovation 
of the molecular theory is mostly due to aggressive ex- 
tension of the scaled-particle-theory concepts. In deed, 
the success of the scaled-particle theory extension d 10 ! 16 ! 
relieves the most serious objections to the PC theory 
which can then be seen as merely an approximate the- 
ory of liquid solutions. Thus applied to molecular-scale 
problems, the PC theory has about the same concep- 
tual status as the distinct minimal scaled-particle theo- 
ries. Of course, the PC theory (and the information the- 
ory approaclP-^ which is a hybrid) addresses hydrophobic 
interactions which have not been addressed by scaled- 
particle approaches. 

Here we address the issue of B 2 from the renovated 
point-of-view. We find that B 2 for atomic-size hard 
spheres in water is attractive (B 2 < 0), and more at- 
tractive with increasing temperature (AB 2 /AT < 0). 



II. THEORY 

We seek the cavity correlation function 

VAA(r) = exp [u A A(r)/k B T}g AA (r) (2) 
for atomic-size hard spheres on the basis of the potential 



1 PC 

distribution theorem (or test particle) approac h 1 18 ! 19 ! 



y AA (i,2) = 



1,2 



(3) 



This formula is cast for evaluation on the basis of trial 
placements of hard-spheres at specific points, here the 
two points (1,2). Of course, we expect ?/aa(1,2) — 
2/aa(?") to depend only on the magnitude r of the dis- 
placement between positions 1 and 2. The notation (. . .) 
indicates the average over the configurations of water 
without the solutes present. AU A 1} = U(N+1)-U(N)- 
U(l) is the binding energy for insertion of an A atom, and 
for the hard sphere case considered here is either zero (no 
overlap with a water oxygen atom) or positive infinity. 

Thus e~ ,9AC/ A ' is an indicator function for permissibil- 
ity of an insertion at a point considered. We will treat 
(2) 

the case that AU\l for two trial placements is additive, 

AU^ + AU^ for the AA atoms considered. 

We rearrange this formula to make the numerical es- 
timation transparent, and this rearrangement overlaps 
a recent discussion of the Kirkwood-Salzburg theory.EO 
Note that the denominator factors of Eq. ^ , being aver- 
ages of indicator functions, are probabilities. Following 
the primitive understanding of conditional probabilities 
p(A\B) — p(AB)/p(B), we use one of those denominator 
probabilities to introduce the expectation conditional on 
permissibility of that initial insertion. Taking the posi- 
tion of that first insertion to be the origin we write 



(4) 



where r is the position of a trial placement relative to 
a permissible insertion. The average indicated in the 
numerator is conditional on the permissible placement 
at 0, though we do not set-up a further notation for 
that. Eq.fi]) expresses the well-known zero-separation 
theorerrPHM 



J/aa(0) = 



1 



(5) 



since the numerator is one (1) under the condition of a 
permissible insertion at 0. 

To estimate the ratio (Eq. Q) we exploit many (n t ) 
trial placements into the system volume V, for each one of 
ti c configurations, c (Fig. [T]). Those trial points will have 
the density n t /V and are statistically uniform. Out of 
the nt trial points, a smaller number n s (c) are permissible 
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FIG. 2. Cavity correlation functions for hard-sphere solutes 
ion liquid water at constant pressure p = 1 atm, and four dif- 
ferent temperatures. The spheres have van der Waals radius 
of 0.17 nm and distance-of-closest-approach to a water oxygen 
atom of 0.31 nm. The dots on the left vertical axis are values 
calculated independently on the basis of the zero-separation 
theorem. 



was 0.31 nm, corresponding approximately to the case of 
an Ar solute. Larger solutes would make the present cal- 
culations prohibitively difficult. The results for In y aa(t) 
(Fig. [2]) operationally satisfy the zero-separation the- 
orem, show strong hydrophobic attraction as short- 
distances, and solvent-separated hydrophobic attraction 
qualitatively consistent with the PC theory. 
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placements, and we estimate the denominator of Eq. Q 
witn E c "s(c)/ (n c n t ) = n s /n t . 

We expect (n t — 1) AV/V of those trial placements to 
land in a volume AV which is a thin shell of radius 
r > surrounding a permissible insertion. Let's denote 
by An s (r;c) the number of permissible placements ob- 
tained in the shell for configuration c. We estimate the 
numerator of Eq. Q as 



J2 C An s (r;c) 
(nt - 1) n c AV/V 



An s (r) 
(n t - 1) AV/V 



Combining these results, we have the estimate 



V 



y A A(r)AV = An s (r) 



(6) 



(7) 



when n t —¥ 00. This is the same formula as if the permis- 
sible insertions were actual particles though that is not 
the case here. 

The formula Eq. Q is operationally consistent with 
the zero-separation theorem Eq. ([5 1 according to the fol- 
lowing argument. Consider a small volume element sur- 
rounding the position r = that is known to be a permis- 
sible placement. We expect that all trial placements in 
this region should be permissible so An s (0) « n t AV/V 
on the right-side of Eq. Q. Therefore, j/aa(0) ~ n t /n s 
which is the operational content of Eq. ([5| . 



III. RESULTS AND DISCUSSION 

Using standard methods (detailed below), this ap- 
proach was implemented for the case that the hard- 
sphere distance-of-closest-approach to an oxygen atom 



FIG. 3. Radial distribution functions for hard-sphere solutes 
in liquid water at constant pressure p — 1 atm, and four dif- 
ferent temperatures. The spheres have van der Waals radius 
of 0.17 nm and distance-of-closest-approach to a water oxy- 
gen atom of 0.31 nm. The prediction of the PC theory^ at 
T = 300 K for the contact value, <?aa(»" = 0.34 nm), is shown 
by the star. The contact values determined here are consistent 
with the information theory modeP^ result published previ- 
ously, and more than twice larger than the predictions of the 
PC theory. These contact values are higher for higher temper- 
atures, indicating stronger hydrophobic contact attractions at 
higher temperatures. 

The radial distribution function gAA( r ) = J/aaM f° r 
r > 2x 0.17 nm = 0.34 nm (and zero otherwise, Fig.|3|, 
scrutinizes these properties more closely. The contact 
values, gAA(r = 0.34 nm), determined here are consis- 
tent with the information theory model result published 
previously,^ and more than twice larger than the pre- 
dictions of the PC theory. The contact values are higher 
for higher temperatures, indicating stronger hydrophobic 
contact attractions at higher temperatures. The contact 
values of the PC theory also increase with increasing 7^ 
but those increases are small, and the PC contact values 
are sufficiently different from the numerical values found 
here that the small increases are not interesting. 

It is commonplace for simulation calculations to deter- 
mine the temperature dependences of the AA potentials 
of mean force uaa(?") = — k^T lngAA(f)- Some exam- 
ples are^2H22l. We note below that the recent work of 
that type supports our conclusions here (Fig. [3j). But 
that helpful recent work also addresses our present prob- 
lems only indirectly because it does not test the statistical 
mechanical theory, nor isolate aspects of interactions of 
different physical type, nor obtain B2 for those purposes. 
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Based on the experience, the measurable osmotic second 
virial coefficients are more subtle, and despite interest- 
ing suggest ions w AA (r)'s have not been measured for 
molecular-scale hydrophobic solutes. 

The B 2 integrals (Fig. [1]) provide the solution thermo- 
dynamic assessment of these distribution functions. The 
suggested B 2 values are decidedly attractive (negative) 
for T = 300 K and become more attractive at higher 
temperatures. The biggest negative contribution is as- 
sociated with contact hydrophobic attractions. Solvent- 
separated hydrophobic attractions near R w 0.7 nm are 
distinct but smaller than contact hydrophobic interac- 
tions. 

The significance of the solvent-separated hydropho- 
bic interactions has been much discus sed f ollowing the 
ground-breaking work of Pangali, et a/! 32 * 33 *. Those sim- 
ulations treated Lennard- Jones (LJ) model solutes some- 
what similar to Kr or Xe solutes. Simulation results 
(with LJ attractions) were compared with PC theory (not 
treating LJ attractions though modified for continuous 
repulsive interactions) . In view of the introductory dis- 
cussion, further analysis of the role of attractive interac- 
tions should be helpful. But additionally, since the simu- 
lation did not determine the additive constant to idaaW, 
the comparison proceeded after matching the two results 
at their minimum value s, a c onvenient choice that has 
been sometimes followed ! 14 * 34 ^ This comparison can give 
the impression that the non-matched solvent-separated 
hydrophobic interactions are unusually variable or sig- 
nificant. If, for the purposes of comparison, the Waa(t) 
were matched at solvent-separated radii, then the Pan- 
gali, et al. results show stronger contact hydrophobic at- 
tractions than does the PC theory, qualitatively in agree- 
ment with the present work. Recent molecular dynamics 
simulations for Xe(aq) or CH^aq) pairs^^ni agree qual- 
itatively with the present results and thus support this 
conclusion. Still the role of longer-ranged attractive in- 
teractions remains to be studied. 

Inclusion of longer-ranged attractive interactions, i.e., 
London dispersion interactions, can change these B 2 val- 
ues and trends depending on the balance of solute-solute 
and solute- water interactions!^ 1 -^ Conclusive information 
for the case that more general interactions are absent, as 
with this work, is crucial to justifying further analysis 
of more general interactions. The LMF theory^ is one 
promising suggestion for how to proceed with inclusion of 
longer-ranged interact ions, has commonalities with ear- 
lier intuitive proposals}^ 1 -^ and deserves further develop- 
ment. 

The temperature dependence observed here has gen- 
erally been considered counter-intuitive, and is some- 
times referred to as an inverse temperature dependence. 
An explanation why this behavior might be considered 
counter-intuitive is the following^ hydrophobic associa- 
tion is typically rationalized as clumping of inert solution 
inclusions due to specific structuring of their hydration 
shells. It might be guessed that the specific structuring 
should be more significant at lower temperatures, so per- 



haps the hydrophobic association should be stronger at 
lower temperatures, perhaps even more important yet in 
super-cooled water. 

Hydration-shell structuring surely is an important fac- 
tor in hydrophobic interactions. What this argument 
does not address is the distinctive equation of state of 
liquid water. Some well recognized peculiarities occur at 
higher-than-physiological temperatures; for example the 
compressibility minimum occurs at 46 C, under these low 
pressure conditions. The eventual statistical mechanical 
explanation^ of the similarly counter-intuitive entropy 
convergence hydrophobic phenomenon (at T w 130 C) 
depended firstly on proper involvement of the actual 
equation of state of liquid waterP^l Indeed the calcula- 
tions here model that specific equation of state also. 
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FIG. 4. Running integral B2[R] (Eq. |TJ) for assessment of 
the R — > 00 value. The dotted curve for R < 0.34 nm is 
the positive contribution from the hard-core of §aa {r) that 
is common to all here. The suggested R — > 00 values are 
decidedly negative (attractive) and become more attractive 
at higher temperatures. The biggest negative contribution 
is associated with contact hydrophobic attractions. Solvent- 
separated hydrophobic attractions near R « 0.7 nm are dis- 
tinct but smaller. The predictions of the PC theory for B2 in 
these circumstances are repulsive (positive). 



IV. CONCLUSIONS 

This work tests the statistical mechanical theory of hy- 
drophobic interactions, isolates consequences of excluded 
volume interactions, and obtains B 2 for those purposes. 
Contact hydrophobic interactions between Ar-size hard- 
spheres in water are significantly more attractive than 
predicted by the Pratt-Chandler theory. The correspond- 
ing B 2 results for atomic-size hard spheres in water are 
attractive (B 2 < 0), and more attractive with increasing 
temperature (AB 2 /AT < 0) in the temperature range 
300 K< T <360 K. This information is essential for fur- 
ther development of the molecular-scale statistical me- 
chanical theory of hydrophobic interactions. 
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V. METHODS 

The GROMAC^H package and the SPC/lP mo d e l 
was used to simulate liquid water. This simulation 
adopted the isothermal-isobaric ensemble at four differ- 
ent temperatures (T= 300 K, 320 K, 340 K and 360 K). 
A Nose-Hoover thermostat maintained the temperature 
and a Parinello-Rahman barostat was used to establish 
the pressure at 1 atm. Bonds involving hydrogen atoms 
were constrained by the LINCS algorithm. Conventional 
periodic boundary conditions and particle mesh Ewald, 
with a real-space cut-off at 1 nm, was used to treat long- 
range interactions. Simulation cells containing 5xl0 3 
randomly placed water molecules were created utilizing 
PACKMOL 40 to match the experimental density approx- 
imately. After lxlO 4 steps of energy minimization and 2 
ns of density equilibration, trajectories of 20 ns (sampled 
1/ps) were obtained at each temperature. Each simula- 
tion frame was analyzed for cavities based on n t — 2 x 10 5 
trial placements with a distance of closest approach to an 
oxygen atom of 0.31 nm. Successful placements can be 
considered as hard spheres of radius 0.17 nm correspond- 
ing approximately an Ar atom. The cavity analysis is 
about an order-of-magnitude more computational effort 
than the generation of the molecular dynamics trajecto- 
ries. 

Considering the integrand of Eq. ([!]) , #aa ( r ) ~ 1 in 
the thermodynamic limit. For a fixed particle numbers 
in a finite system, that subtracted uncorrelated feature 
is less than 1 and the correction is 1 , " . vanishing 

in the thermodynamic limit. In the present work, we do 
not have fixed numbers of A particles. 
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